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Abstract 

Tau neutrinos interacting inside the Earth produce r leptons which thereafter 
can decay inside the atmosphere. The propagation of extremely energetic v T J s and 
r's through the Earth is studied by means of a detailed Monte Carlo simulation, 
taking into account all major mechanisms of v T interactions and r energy loss as 
well as decay modes. The rates of r's emerging from the Earth are determined as a 
function of r's energy for several cosmic neutrino models. 
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1 Introduction 



Cosmic neutrinos give a unique opportunity to open a new observational window in the 
field of Astrophysics and Cosmology. The detection of such particles, not deflected or ab- 
sorbed during their intergalactic path, could be revealing for the identification of sources 
of Ultra and Extreme High Energy Cosmic Rays (UHECR and EHECR). The detection of 
such ultra and extreme high energy neutrinos is a challenge as it demands a very large sen- 
sitive area. The experiments of the upcoming generation plan to use underwater-underice 
Cherenkov detectors |lj and large field-of-view atmospheric fluorescence detectors ||. 
One of the best evidence of cosmic neutrinos would be the detection of upstream show- 
ers/particles emerging from the Earth. For this kind of events the atmospheric muon 
and primary charged cosmic ray background would be completely suppressed. On the 
other hand this signature can hardly be observed at extreme energy because the rise of 
weak cross sections entails the opacity of the Earth with respect to neutrino propaga- 
tionlli. 

The problem of muon neutrino propagation through the Earth can be solved by taking 
properly into account neutrino-nucleon scattering inside the Earth, either by Monte-Carlo 
simulation || or by approximate iterative solution of the relevant transport equation PJ . 
It has been recently pointed out in references f| @ |K| 01 that the behaviour of r- 
neutrinos, whose existence should be guaranteed in a neutrino-oscillation scenario, should 
be significantly different from and v e . Whilst muon and electron neutrinos are prac- 
tically absorbed after one charged current (CC) interaction, the r lepton created by the 
v T CC interaction may decay in flight before losing too much energy, thereby generating 
a new v T with comparable energy. Hence ultra high energy r-neutrinos and r's should 
emerge from the Earth instead of being absorbed. For a correct evaluation of the energy 
spectrum of u^s and r's emerging from the Earth, one has to properly take into account 
v T interactions as well as r energy loss and decay. An analytical approach similar to that 
used in f|] for z/ M has been proposed by S. Iyer et al. [|l(| neglecting r energy loss. This 
approach holds as long as energy does not exceed 10 16 eV because above this energy r 
interaction length becomes comparable with r decay length ||. 

In reference || a detailed Monte Carlo calculation of v T — r system propagation through 
the Earth has been performed for energy up to 10 20 eV including the r energy loss contri- 
bution, with special emphasis on the initial v T spectrum deformation as a function of the 



zenith angle of the emerging particle. As pointed out in ref. [[IT]] |T2[ [|T3| the tau leptons, 
created by CC v T interactions inside the Earth, because of their relative long decay length 
at high energy, could emerge from the Earth surface and eventually decay inside the at- 
mosphere (fig. HI). Such kind of events could be detected by atmospheric shower detectors 
as upward going showers. In reference JTTJ an overview of the physics processes involved 
in case of r's emerging from Earth and mountains is presented, suggesting the possibility 
of the detection of such events. In reference [|EJ the rate of expected events in the Auger 
detector is given on the basis of a Monte Carlo simulation while in reference | I3fl , and 
more extensively in the present work, the results of a detailed Monte Carlo simulation are 
presented in a general framework, giving the rate of emerging r's for unit Earth surface 
above a given minimum energy, from 10 14 eV to 10 22 eV and for a large variety of neutrino 
fluxes. The results are presented considering two different calculations of r photonuclear 
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cross section and comparing directly the rate of emerging r events to the rate of Extensive 
Air Showers (EAS) produced by neutrinos interacting inside the atmosphere. 



2 Lepton interactions 

Deep inelastic neutrino-nucleon scattering is the dominant interaction of energetic neutri- 
nos into conventional matter. Charged and neutral current differential cross sections |J 
used for this simulation |J have been calculated in the framework of QCD improved par- 
ton model. Since a large contribution to cross sections comes from the very low x region, 
laying beyond present accelerator domain, an extrapolation of parton distribution at very 



low x (x ~ 10 ) is necessary. We have used the parton distribution set CTEQ3-DIS fl4 



available in the program library PDFLIB [15], with NLO-DGLAP formalism used for the 



evolution with the invariant momentum transfer between the incident neutrino and out- 
going lepton (— Q 2 ) and assuming for very low x the same functional form measured at 
x = O(10~ 5 ). Even if more sophisticated approaches for low x extrapolation have been 
developed using dynamical QCD [|1(J the results for cross section calculations (fig. 0) do 



not differ more than 10% from the approach taken here with CTEQ3-DIS plus "brute 
force" extrapolation [|] |L6 . 



The electromagnetic interactions of muons with matter, at the considered energies, are 



dominated by radiative processes rather than ionisation |I7| . Cross sections for elec- 
tromagnetic radiative processes of r are lower than muon's, yet radiative interactions 
remain the dominant process for r energy loss. The cross sections used for radiative 



interactions of r leptons are based on QED calculation for bremsstrahlung [18[ and for 



direct pair production [Hi] while for photonuclear (PH) interactions references 



have been used. For all aforementioned processes we have implemented stochastic inter- 
actions for v = (Ei — Ef)/Ei > 10~ 3 and a continuous energy loss approximation for 
v = (Ei — Ej)/Ei < 10~ 3 , where Ei and Ef are r energies before and after the interaction 
respectively. 

It is worth mentioning that bremsstrahlung cross section scales as the inverse square of 
lepton mass m; whereas direct pair production and photonuclear interaction cross sections 
approximately scale according to 1/m/ p2"| . As a consequence the dominant processes 
for r lepton energy loss are direct pair production and even more photonuclear interac- 
tion rather than bremsstrahlung photon radiation [f2lfl . The calculation of cross sections 
in case of PH interactions has been explicitely performed by Bezrukov and Bugaev [2C 
only for muons while the recent calculation performed by Iyer et al. |yj includes the 
r lepton. The two calculations, as far as r lepton is concerned, differ by almost one 
order of magnitude at extreme energy (fig. ^) and therefore the rate of emerging r's is 
expected to be sizely affected. The r decay has been simulated by using the TAUOLA 
package [p3fl . The Monte Carlo simulation has been performed following neutrinos and 
charged leptons along their path inside the Earth. The Earth model considered is the 
preliminary Earth model of ref. [p4[ , applied to a perfectly smooth Earth surface with 



the outer layer made either of "standard rock" or water. For what concerns the Earth 
composition we have used "standard rock" [Z = 11, A = 22) for the crust and the mantle 
and iron [Z = 26, A = 55.8) for the core. 
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3 Neutrino induced r's inside the Earth 



3.1 Effective depth 

Tau neutrinos interacting in proximity of Earth's surface can produce r's which thereafter 
may survive until they get into the atmosphere. The production of r's which emerge 
from the Earth is driven by two, not completely separated, mechanisms. The first is 
the propagation of neutrinos through the Earth, which determines the amount and the 
energies of neutrinos approaching the crust just below the Earth's surface. The second 
is the production and propagation of r's in the crust just below the Earth's surface. In 
the next section we will show the results of the simulation concerning the whole "story" 
of propagation of leptons, while in this section we focus the attention only on the second 
mechanism. Here we want to evaluate the thickness of the Earth's Crust to be considered 
active for production of emerging r's as a function of the energy of the r neutrino as 
it approaches the emerging surface. Hence no propagation through the whole Earth is 
considered in this section. The results given here can be useful also for calculations of the 
rate of horizontal r's emerging from mountains. 

As a first approximation, the r decay length 7/3cr can be used to estimate the effective 
depth, defined as the thickness of Earth's Crust to be considered as a source for such 
emerging r's. 

An accurate calculation of the effective depth has to properly take into account the 
differential cross sections describing the neutrino interactions and tau energy loss and 
decay. The total neutrino cross section influences the number of produced r's but not 
the effective depth. We have simulated a number N\ nt of v T CC interactions induced by 
vertical neutrinos and uniformly distributed along a column of Earth's Crust within a 
depth L. One can then define the effective depth as the limit: 

L eff = r lim Le{L) 

where e(L) is the number of emerging r's divided by the number of simulated v interactions 
N" nt . The effective depth is shown in fig. § where it is compared with the r decay length 
with E T = E v . At low energy, L e jf/E v slowly increases with E v because the fraction of 
energy carried by the r increases with energy in the CC v T interactions. Above ~ 10 17 eV 
the effect of r energy loss becomes relevant and the energy degradation of r's prevents 
a sizeable fraction of them from reaching the surface. Consequently, L e ff/E u drops and 
L e ff roughly approaches the value of r radiation length in the outer layer of Earth's Crust 
10 km for the highest energies and PH in ref. |20f ). 



3.2 Effective aperture 

In this section the calculation of the fluxes of r's emerging from the Earth has been 
performed by means of a full simulation of v T — r system propagation through the Earth. 

Fig. H shows the result of the simulation for an isotropic flux of vfs impinging on the 
Earth with an energy E v = 10 20 eV. Most r's with zenith angles larger than 94° emerge 
with energies 10 16 eV-10 17 eV for which, assuming PH interaction in ref. [20|, r decay length 
becames comparable with its radiation length. 
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For zenith angles corresponding to a path length in the Earth larger than v T interaction 
length in the Crust [L\ nt ~ 70 km for E u = 10 20 eV) , neutrinos interact far away from the 
Earth surface, and the produced r's lose energy and decay before getting out. Owing to 
this mechanism, the energy of the v T -r system decreases as zenith angle increases (fig. |5]) 
and the probability for a v T to emerge as a r drops at large zenith angles. For v T energy 
lower than E v = 10 20 eV the behaviour of the u T -r system is similar but shifted to larger 
values of zenith angles. 

The number of emerging r's with energy larger than E min (minimum energy) from an 
Earth infinitesimal surface dS can be calculated in a unidimensional approach as : 



dN T (E T > E„ 



d 2 N 
dE„ dSl 



(E„) P Ut ^ t (cos(G), E v , E min ) dS \cos(Q)\ dQ dE, 



where d ^ N dn (E u ) is the cosmic r neutrino flux, B is the zenith angle of the emerging 
particles, P^ T (cos(9), E u , E min ) is the probability that a u T impinging on the Earth 
surface will emerge as a r with an energy E T > E min . Assuming a spherical shaped Earth 
and an isotropic neutrino flux, Eq. [I] turns into : 

d 2 N 



N T (E T > E n 



S 



dE„ dQ 



(E v ) J P Ut ^ t (cos(Q), E u , E min ) \cos{Q)\ dQ dE v 



(2) 



where S is a portion of Earth surface. 

One can now define an effective aperture A e ff(sr) containing all the physics of leptons 
propagation : 

A eff (E u ,E min ) = / P VT _ T (cos{Q),E u ,E min )\cos{Q)\dQ 



so that Eq. reads : 



N T (E T > E n 



S 



Ei/^Ejy 



d 2 N 
dE„ dVl 



(E u ) A e ff(E l/ , E min ) dE u 



(3) 



The effective aperture A e ff does not depend on neutrino fluxes or detector performance 
and allows a simple and straightforward calculation of the number of emerging r's for 
any cosmic neutrino flux and detector sensitive area. It must be pointed out that in the 
definition of P Ut ^ t (cos(Q), E u , E min ) we also constrain the emerging r to decay within an 
atmospheric slant depth larger than 1000 g/cm 2 . This condition guarantees the full size 
formation of an EAS after r decay. 

The results for the effective aperture by using Bezrukov and Bugaev PH cross section 
are shown in figs. ^ and [7|. They seem to be in fairly good agreement with those in 



ref. [0] for the Auger detector, though the contribution of the effective Auger sensitivity 
can not be completely unfolded from quoted results using available informations. 

An effective aperture A u e e ^ for downwards going v e CC interactions in the atmosphere 
is also shown in figs. |6] and ^. For this kind of events, the energy of produced EAS is 
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the same as u e . For an isotropic neutrino flux, the number of EAS induced by u e CC 
interactions within the atmospheric mass overhanging a portion of Earth surface S is 
given by : 

f d 2 N 

Nea S (E eas > E min ) = S / -——(E v )acc(Ev)N AP 2TrdE v 

JE v >E mm dE u d\l 

where p is the atmospheric pressure (column density) at ground level, N A is the Avogadro 
number, d ^ N dn (E v ) is the cosmic electron neutrino flux and a C c{E v ) is the neutrino CC 
cross section on isoscalar target. Again, we can define an effective aperture : 

A u e } f {E u ) = 2nN A pa cc {E u ) 
such that the number of u e CC interactions can be written as: 

N% AS {E EAS > E min ) =S f ^L(^) A ?fA E ») dE » 

A y e e ^{E v ) increases with energy following the rise of neutrino cross section. At low 
energies (E u <C 10 18 eV) the rise with energy of A e ff for emerging r's is steeper than the 
rise of A v e £ ^{E v ) because both the increase of v cross section and the increase of r decay 
length contribute to enlarge P u ^ T (cos(Q), E u , E min ). A e ff{E v ) keeps increasing as long 
as the increase of v cross section enhances the production of r's in the Earth. Around 
E p ~ 10 18 eV the cross section of neutrinos is so high that the probability for neutrinos 
to have at least one interaction inside the Earth is close to one for most zenith angles 
and the point of first v T interaction in the Earth is likely to be far from the emerging 
point (this latter effect is responsible for the slow decrease of A e ff between 10 19 eV and 
10 20 eV in fig|J). Moreover, the increase of r decay length with energy prevents some of 
the extreme high energy events from producing a detectable shower into the atmosphere 
(only r's decaying within 1000 g/cm 2 slant depth are retained). Due to these mechanisms, 
the value of A e ff at extreme high energies shows an almost flat dependence on energy 
and strongly depends on the required E min . 

The effect of decay length selection at extreme energies is shown in fig [| while fig. [F] 
shows the A e ff in case of an outer layer of Earth's Crust made of water and 3 km thick. 
Due to the reduction of target mass, the values of A e ff in fig. ^ are lower, at low energies, 
than those in fig. |6|. For energies E v > 10 18 eV, where A e ff mainly depends on r radiation 
length in the outer Earth's Crust, the values in fig. [7] are larger than those in fig.|6]. 

The angular distributions of emerging r's are shown in fig. || As the initial neutrino 
energy increases, the maximum of the angular distribution moves towards 90° zenith 
angle. 



In fig. 10, the effective apertures calculated using the recent evaluation of r PH cross 



section given in ref. |2T[] are compared with the effective aperture of fig. |6]. The higher 



value of r energy loss (fig ||) gives rise to a suppression of emerging r's for the highest 
energies. Moreover, A e ff decreases at extreme neutrino energy because of the increase 
with energy of r PH radiation length according to fig. |[ 
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3.3 Fluxes of emerging r's 



According to the hypothesis of — u T neutrino oscillations suggested by the Super- 
Kamiokande experiment |25], half of the original cosmic i^'s appear as z/ T 's after their 
intergalactic path towards the Earth |26|. 

Within this scenario, we have used the same fluxes for parent v T and v e : = <& Ue = 
^initial 1 2 anc [ we h ave calculated the rate of emerging r's by using equation (|3]). 

Fig. [IT], [ITJ show the flux of emerging r's compared with the rate of v e CC interactions 
in the atmosphere overhanging the unitary surface. The signal induced by the decay 
of emerging r's dominate for E m i n ^5x 10 17 eV while for higher minimum energy E m i n , 
Earth matter strongly suppresses the upgoing r fluxes. It must be emphasized that, in this 
analysis, we have not taken into account the difference between induced EAS generated 
by different r decay channels. It is worth pointing out that leptonic decays r~ — > n~ v n v T 
do not produce any EAS. 



4 Conclusions 

A detailed study of v T -r leptons propagation through the Earth indicates the likely exis- 
tence of a quite intense flux of upgoing r's emerging from its surface. The rate of EAS 
induced by r's decay is larger than the rate of z/ e 's induced EAS for minimum energy 
E m in < 5 x 10 17 eV, while it drops for larger values of minimum energy. The rate of 
emerging r's at extreme energies strongly depends on photonuclear contribution to en- 
ergy loss, whose calculation is still controversial as it requires non-trivial extrapolations. 
The detection of such r induced events is challenging at a detector's energy threshold 
E t hr ~ 10 19 eV even for new generation experiments, while for E t hr «C 10 18 eV the de- 
tection seems to be realistic for detectors capable of recording signals induced by very 
inclined showers with sensitive areas 3> 10 3 km 2 . 
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Figure 1: Schematic drawing of atmospheric showers induced by neutrinos. Upward going 
IV 's could be detected as upward going showers induced by r decay. 
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Figure 2: Total v T CC and NC cross sections on isoscalar nucleon target as a function 
of energy. Continuos line shows the neutrino CC cross section; dotted line shows the 
antineutrino CC cross section; dashed line shows the neutrino NC cross section and dotted- 
dashed line shows the antineutrino NC cross section. CTEQ3-DIS parton distribution has 
been used for the calculation. 
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Figure 3: (3 — — dE ^ x I in case of standard rock (A=22, Z=ll) for r PH interactions by 



S. Iyer et a/. |2T[ and as it comes out using formulae of double differential cross section 
of Bezrukov and Bugaev [2C] with r mass. 
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Figure 4: Effective depth for generation of emerging r's (PH cross section from ref. 
and outer Earth layer made of standard rock). The dotted line shows the expected result 
for no interacting r's created with the same neutrino energy and decaying after one r 
decay length. 
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Figure 5: Scatter plot of final energies of emerging r's {Ef) versus emerging zenith angles. 
The simulation has been performed at E v = 10 20 eV for an isotropic flux. PH cross section 



is from ref [20| and the standard rock is considered for the outer Earth layer. 
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Figure 6: Effective aperture for different minimum energy E r 
sidered for the outer Earth layer. PH cross section is from ref pi] 



and standard rock con- 
Dotted line shows 



the effective aperture for events induced by downward going v e CC interactions in the 
atmosphere. 



15 



oo 10 



.a- 



< 



Water 



-4 



-5 
=- 



-6 
=- 




14 15 1i 

10 10 



• 


E rnin 


= 1*' 


9 eV 


* Emir, = 8* ' 


,7 eV 





^min 


= 8*' 


13 eV 


o E m!n = 8*1 


17 eV 




^min 


= 8*' 


14 eV 


■ E ■ =3*1 


0' 8 eV 




Ernin 


= 8*' 


15 eV 


a E min = 8*1 


0' 8 eV 


T 


E m ;„ 


= 8*' 


O 16 eV 


★ E ■ =3*1 


,9 eV 




I 




i 


□ v e atm. 










1 7 




1 9 




22 




io ,s 


10 


io 20 io 21 


10 



E„ e 



Figure 7: Effective aperture for different minimum energy E min and 3 km water considered 



for the outer Earth layer. PH cross section is from ref ||20|| . Dotted line shows the effective 
aperture for events induced by downward going v e CC interactions in the atmosphere. 
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Figure 8: Effective aperture for minimum energy -E m j n = 3 x 10 19 eV, standard rock 
considered for the outer Earth layer and PH cross section from ref |j20 |. The effect of 
different selections on the slant depth of the decay point is shown. 
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Figure 9: Angular distribution of emerging r's for several primary neutrino energies and 
several minimum energy E min . PH cross section is from ref |2(J and standard rock is 
considered for the outer Earth layer. 
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Figure 10: Effective aperture for different minimum energy E min and standard rock con- 
sidered for the outer Earth layer. Calculations using Bezrukov and Bugaev r PH cross 
section |20] (curves labeled as PHI) are compared with calculations using S. Iyer at al. 
cross section |H| (curves labeled as PH2). Dotted line shows the effective aperture for 
events induced by downward going v e CC interactions in the atmosphere. 
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Figure 11: Rates of emerging r's versus minimum energy : Bezrukov-Bugaev PH cross 
section [2(| and standard rock considered for the outer Earth layer (continuos lines), 



Bezrukov-Bugaev PH cross section [2D| and 3 km water for the outer Earth layer (dotted 
lines), S. Iyer et al. PH cross section [21| and standard rock for the outer Earth layer 
(dashed lines). Dashed-dotted lines show the rate of v e CC interactions in the atmosphere 
above the unit Earth surface. Each plot refers to a different Bottom- Up model of cosmic 
neutrino flux: (a) AGN (A) from Mannheim [27]; (b) AGN (B) from Mannheim [27| ; 
(c) GRB from Waxman e Bahcall |58] and Vietri extension P9| 
Protheroe and Johnson 130 
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Figure 12: Rates of emerging r's versus minimum energy : Bezrukov-Bugaev PH cross 
section [2(| and standard rock considered for the outer Earth layer (continuos lines), 



Bezrukov-Bugaev PH cross section [2~U| and 3 km water for the outer Earth layer (dotted 
lines), S. Iyer et al. PH cross section [21| and standard rock for the outer Earth layer 
(dashed lines). Dashed-dotted lines show the rate of v e CC interactions in the atmosphere 
above the unit Earth surface. Each plot refers to a different Top-Down model of cosmic 
neutrino flux: (a) Topological Defect from Protheroe and Stanev (b) Topological 
Defect from Sigl et al. |52]] ; (c) Topological Defect from Sigl (d) Supermassive Relic 
Particle from Kalashev et al. 
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